^ 


Richard  A.  Reinhardl 

OPERATIONAL  ANALOGUES 
FOR  KINETIC  STUDIES. 


TA7 
.U6 
no. 15 


library 

U.  S  Postgraduate  School 


UNITED  STATES 
NAVAL  POSTGRADUATE  SCHOOL 


OPERATIONAL  ANALOGUES  FOR  KINETIC  STUDIES 

by 

RICHARD  A.  REINHARDT 
Assistant  Professor  of  Chemistry 

GILBERT  F.  KINNEY 
Professor  of  Chemical  Engineering 


RESEARCH  PAPER  NO.  15 
November  1957 


«s  ^ 


<p 


0 


& 


OPERATIONAL  ANALOGUES  FOR  KINETIC  STUDIES 


by 

RICHARD  A.  REINHARDT 

Assistant  Professor  of  Chemistry 

GILBERT  F.  KINNEY 
Professor  of  Chemical  Engineering 


Research  Paper  No.  15 
UNITED  STATES  NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 
November  1957 


US.  NAVAL  POSTGRADUATE  SCHOOT 


FOREWORD 

This  paper  was  presented,  upon 
invitation  of  the  symposium  chairman, 
as  part  of  the  Beckmann  Award  Sym- 
posium at  the  131st  meeting  of  the 
American  Chemical  Society  on  8   April 
1957  in  Miami,  Florida. 


ABSTRACT 

Programming  of  the  operational 
analogue  computer  for  the  rapid 
solution  of  problems  dealing  with 
complex  chemical  reaction-rate 
mechanisms  is  described.  As  an 
illustration,  an  enzyme-catalyzed 
reaction  is  considered  and  it  is 
shown  that  computer  results  agree 
well  with  experiment,  but  not  with 
the  steady- state  approximation 
customarily  made  for  this  system,. 
It  is  pointed  out  that  the  oper- 
ational analogue  computer  is  a 
relatively  simple  and  inexpensive 
set-up,  and  that  programming  prin- 
ciples, as  developed  and  described 
here,  make  it  applicable  to  a  wide 
variety  of  theoretical  and  industrial 
problems. 
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An  operational  analogue  computer  is  an  assemblage  of 
independent  high-gain  DC  amplifiers,  with  accessory  resistors, 
capacitors,  and  potentiometers,  and  a  peg-board  arrangement 
for  making  connections.   Each  time-dependent  variable  is 
represented  on  the  computer  as  a  voltage,  and  the  amplifiers 
can  either  add  together  several  such  voltages,  or  else  add 
and  integrate.   It  is  the  purpose  of  this  paper  to  show  how 
an  analogue  computer  of  this  type  can  be  used  for  investigating 
the  mechanism  of  a  complex  chemical  reaction.   The  unravelling 
of  such  a  mechanism  would  then,  for  example,  permit  ready 
scale-up  from  laboratory  to  plant  size  equipment,  and  thus  be 
of  valuable  assistance  in  the  selection  of  optimum  size  reactors, 

In  many  actual  applications  of  chemical  kinetics,  the  rate 
law  corresponding  to  a  given  mechanism  is  approximated  by  the 
use  of  simplifications  suitable  to  the  problem  at  hand  -  such 

as  the  steady-state  hypothesis  customarily  required  in  a  large 

3 

number  of  chemical  mechanism  studies.     When,  however,  the 

concentration. of  each  species  (including  intermediates)  is 
desired  as  a  function  of  time,  especially  in  the  early  stages 
of  reaction  before  the  steady  state  is  established,  it  becomes 
necessary  to  use  the  complete  integrated  form  of  the  rate  law. 
This  form  is  also  required  if  the  rates  of  several  steps  are 
of  comparable  magnitude,  so  that  the  concept  of  a  rate- 
determining  step  is  not  significant.   These  complete  expressions 
are  extraordinarily  complex  and,  indeed,  in  many  cases  can  be 
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evaluated  only  by  numerical  methods »   This  is  especially  true 
for  reactions  involving  bimolecular  or  trimolecular  steps  in 
the  mechanism,, 

An  alternative  treatment  for  the  complete  rate  law  is  that 

of  working  with  the  set  of  simultaneous  differential  equations 

6" 

directly  obtainable  from  the  mechanism,,   Each  of  these 

equations  is  of  first  order  and  first  degree,  though  non- 
linear for  bimolecular  or  trimolecular  steps.   (It  should  be 
noted  that  even  though  the  reaction  is  of  higher  than  first 
order  in  the  chemical  sense,  each  differential  equation  is 
mathematically  of  the  first  order,) 

Consider,  for  example,  a  reaction  of  two  consecutive 
unimolecular  steps  J 

A_J_^B l^C  (1) 

Selecting  the  rate  equation  for  component  B  as  an  example, 
and  using  lower-case  letters  to  represent  concentration, 

db   =  k,a  -  k0b  (2) 

dt       L  * 

Defining  the  operator  l/P  by  the  equation 

t 


l/P  [F(t)J  =    f   F(t)dt, 


(3) 


equation  2  may  be  integrated  to 


t 


b  =    r       (k.a  -  k  b)dt  +  bn       (4) 


J 


0 


FIG.  I    THE  INTERGRATING    AMPLIFIER:  (upper) 
ELECTRICAL   CIRCUIT;  (lower) SCHEMATIC. 


or 


-1/P  (k2b  -  k^a)  +  bQ  (5) 


in  which  subscript  zero  refers  to  the  concentration  at 
zero  time. 

The  upper  portion  of  Fig.  1  shows  the  electrical  circuit 
for  an  amplifier  arranged  to  perform  the  operation  indicated 

by  equation  5;  this  is  shown  schematically  in  the  lower 

7 
portion  of  the  figure,  using  conventional  computer  symbols. 

The  amplifier,  with  amplification  -A,  is  shown  as  a  triangle, 

the  grid  input  being  to  the  base  of  the  triangle,  the  plate 

output  from  the  apex.   The  potentiometers,  set  at  a-, and  a^, 

control  the  fraction  of  the  input  voltages,  e^and  e2 (referred 

to  ground)  fed  through  the  input  resistors,  R,and  R?,  to  the 

grid  of  the  amplifier.   For  feed-back  capacitance  C  and  with 

output  voltage  e,  (also  relative  to  ground)  it  follows  from 

Kirchhofffs  Law,  applied  to  the  grid  input,  that 

a  e  -  e      a^e«-  e     „  d(e  -  e)    .    rt      ... 
1  1   g  +   2  2   g  +  C  £.  +  ig  =  0,     (6) 

,    Rl  R2  dt 

where  subscript  g  refers  to  the  grid.   Now  the  grid  voltage 

and  current  are  very  small  for  a  well-designed  high-gain 
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amplifier,  and  from  this  consideration  it  develops*    that, 

in  terms  of  the  operator  1/P,  the  output  voltage  becomes 

6  ■  -1/p(^-  ei  +  ^T  e2)  +  E°'    •  (7) 

where  E  is  the  initial  value  of  the  output  voltage 
(disconnected  when  computer  operation  begins). 

-  k   - 


■ 


e,  _ 


'€> 


e2 @ R2 


FIG.  2     THE   SUMMING  AMPLIFIER  (upper) 

SET-UP   FOR  THE  ENZYME  PROBLEM,  (lower) 


The  similarity  in  form  between  equations  5  and  7  suggests 
the  use  of  such  amplifiers  in  solving  the  type  of  problem 
exemplified  by  equation  1,  The  requirement  is  that  the  vol- 
tages suitably  represent  concentrations  and  the  resistor- 
capacitor-potentiometer  combinations  represent  the  rate 
< 

constants. 

In  many  problems  it  is  also  necessary  to  add  together 
variableso   For  this  purpose  can  be  used  the  summing  am- 
plifier, as  shown  schematically  in  the  upper  portion  of 
Figure  2.   The  arrangement  is  similar  to  that  for  the  inte- 
grating amplifier,  save  that  the  feed-back  is  resistive,  of 
value  R„.  Making  assumptions  as  previously,  the  equation 
for  the  output  voltage  is 

a  a 

e  =  -Rf(  1  e    +   2   e2)  (S) 

Rl         R2 
For  single  input  and  equal  resistors,  the  summing  amplifier 
functions  as  a  sign  changer,, 

If  the  mechanism  under  consideration  contains  steps  of 
molecularity  greater  than  one,  it  is  also  necessary  to  use 
a  function  multiplier,  which  is  an  electronic  device  to 
generate  the  product  of  two  variable  voltages  <>   For  commer- 
cially available  multipliers  the  output,  as  a  function  of 
input  voltages  e-,and  e2,  is  given  by 

e  =  -  me,e2,  (9) 
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where  m  is  a  constant  of  proportionality  (equal  to  0„01  for 
the  equipment  used  in  this  laboratory) »  On  the  schematic 
diagrams  the  multiplier  is  shown  as  a  rectangle  with  an  X 
at  its  center.   Another  auxiliary  piece  of  equipment  is  a 
function  generator,  whose  output  conforms  to  some  arbitrarily 
given  value  for  any  given  input »   This  device  is  useful  for 
non-isothermal  reactions,  where  rate  constants  may  change* 

To  illustrate  the  use  of  computer  technique,  the  following 
mechanism  has  been  selected: 


A  +  B 


k 
C  3 


B  +  D 


(10) 


in  which  B  is  a  catalyst  of  some  sort.  Equation  10  could, 
for  example,  represent  the  Michaelis-Menten  formulation  for 
enzyme  activity,    B  symbolizing  the  enzyme  and  A  the  sub- 
strate.  Assuming  only  A  and  B  present  initially,  the 
operator  equations  for  A,  C,  and  D  become 


a  =  -  1/P 

-c  =  -  1/P 
d  =  -  1/P 


kiab  +  k2^ 


-o)]  ♦  ao 


kiab  +  k2*~°*   +  k3^"cJ 


y-o] 


(ill 

(12) 
(13) 


The  equations  are  written  in  terms  of  -c  to  avoid  the  use 
of  additional  amplifiers  for  sign  changing.   The  equation 
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for  B  can  be  deduced  from  the  stoichiometry  of  the  system: 


■b  = 


bQ  ♦  (-C) 


(14) 


-b  being  used  here  since  the  function  multiplier  changes  sign. 
These  four  equations  suggest  the  use  of  at  least  four 
amplifiers  (three  integrating  and  one  summing)  plus  a  function 
multiplier.   The  arrangement  used  to  solve  this  problem  is 
shown  schematically  in  the  lower  portion  of  Figure  20 

The  specific  problem  chosen  for  illustration  was  the 
enzyme-catalyzed  oxidation  of  an  appropriate  oxygen  acceptor 
by  hydrogen  peroxide,  a  system  which  has  been  studied  by 


Chance 


Typical  values  of  the  parameters  for  this  problem 


are: 
a 


-6  -6  7 

0  =  4  x  10    M,  bQ  =  1x10   M_>ki  =  1x10  ,  k2=  0.2, K  =5 


(t  in  seconds  for  the  rate  constants).,   To  illustrate  the  scaling 
it  will  suffice  to  consider  only  the  amplifier  generating 
a,  the  voltage  representing  a«Let  ^A  =  a/a,<^B  -  b/b 
\   c  =  c/c,  where  each  ^\  is  a  magnitude-scaling  factor,  and 
"a",  b",  and  c  the  voltages  representing  the  respective  con- 
centrations o  Then  equation  11  may  be  rewritten 


a  = 


1/P 


^ 


l0 


(15) 


lr  X  ab  +  k2  X—  (~c) 

To  take  care  of  time-scaling,  a  factor  p  is  introduced  with 
each  feed-back  capacitance,  such  that  p  seconds  on  the 
computer  corresponds  to  one  second  of  actual  time..   Then  by 
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applying  the  electrical  equations  7  and  9  to  this  amplifier, 


=  -i/p 


jJ:  is  +  Jii;   (-c) 

+  a 

RXC                       R2C 

0 

(16) 


where  subscript  1  refers  to  the  portion  of  the  circuit  into 
which  maT  is  fed,  subscript  2  to  that  for  -c0  By  equating 
coefficients  in  equations  15  and  16, 


i  \  -  & 


(17) 


r-lC 


k2  ^C  B  a2  ft  (IS) 

^A    r2c 

These  two  equations^  combined  with  similar  equations  for 
the  other  amplifiers ,  permit  the  proper  selection  of  the 
resistors,  capacitors,  and  potentiometer  settingSo  In  Table 
1  are  shown  values  of  these  components  corresponding  to  one 
case  considered ,  and  in  Figure  3  are  shown  simultaneous 
traces  (from  a  four-channel  recorder)  of  each  of  the  concen- 
trations vs_  time  for  this  case.   Concordance  with  the  results 

1 

of  Chance     are  quite  satisfactory?  e_0go ,  the  maximum 

concentration  of  the  intermediate  C  and  the  time  required 
to  attain  this  maximum  agree  quite  well  with  the  results 
reported  by  himc 

The  assumption  that  there  is  established  a  steady  state 

4, 
in  b  and  c,  after  an  initial  transient  period  leads  to 

the  relationship 


Table  I.   Typical  Values  of  Components  for  the  Enzyme  Problem 


Amplifier  (Species) 


Magnitude- scaling  factor, 
^ ,  x  10* 

Feed-back 

Input  Resistors,  megohms 

Initial  voltage 

Initial  concentration, 
micromoles/liter 


A 

B 

c 

D 

10 

1 

1 

10 

1/4 

•  2  meg 

y 

V* 

1,    10 

O   <~                       O  -L    j 

i,  i 

1 

3* 

95 

0 

0 

3»3 


.95 


Potentiometer  Setting 
Rate  Constant  (t  in  sec*) 


Potentiometer  (Rate-constant) 


1 

■   i — 
2 

3 

loO 

0,02 

0.5 

7 

1x10 

0<,2 

5 

Time-scaling  factor,  p  ,   =10 
Multiplier  scaling  factor,  m,  =  0o01 
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FIG.3     KINETICS   OF  THE  ENZYME  PROBLEM 
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FIG.  4  EFFECT  OF  CATALYST  CONCENTRATION 
ON  STEADY— STATE    ASSUMPTION. 
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k     k 
c  =    2  !  3  ab,  (19) 

kl 

which  applies  strictly  only  when  db/dt  =  dc/dt  =  0„   In 
Figure  4  is  shown  a  portion  of  the  plot  of  c  vs„  ab,  for  a 
=  4  x  10    M,  and  various  values  of  b  ,  as  made  directly 
from  the  computer  on  an  X-Y  recorder „   (note  that  time 
increases  from  right  to  left  in  this  figure.)   The  straight 
line  is  that  computed  from  equation  19»   It  can  be  seen  that 
if  b  is  comparable  with  a  ,  the  steady- state  approximation 
is  not  valid,  even  near  the  completion  of  reaction,.   Only  at 
small  values  of  b  /an  do  the  actual  curves  approach  equation 
19,  although  the  straight  line  could  serve  as  a  rough  average 
for  the  latter  stages  of  reaction,.  Whereas  deductions  such 
as  these  are  readily  made  from  computer  traces,  it  is  indeed 
quite  difficult,  if.  not  impossible,  to  produce  them  by 
ordinary  mathematical  methods* 

For  a  general  kinetics  problem  at  least  one  amplifier 
will  be  required  for  each  component  participating  in  the 
mechanism,  plus  such  others  as  may  be  required  for  sign- 
changing;  one  potentiometer  for  each  rate-constant;  and  one 
function  multiplier  for  each  bimolecular  step  (or  two  for  a 
trimolecular  step,  etc.).   The  speed,  simplicity,  and 
adaptability  of  the  operational  analogue  computer  makes  it 
also  a  useful  device  for  the  determination,  from  observed 
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kinetics  data,  of  the  rate  constants  for  assumed  mechanisms. 
Trial  values  of  the  various  rate  constants  are  used,  by- 
changing  potentiometer  settings  and,  if  necessary,  resistors, 
until  a  fit  with  the  data  is  obtained.   Since  changing 
parameters  ordinarily  takes  only  a  few  seconds,  such  a 
problem  can  be  solved  quite  rapidly. 

It  follows  that  the  independent  variable  of  the  problem 
need  not  be  time,  but  could  also  be  any  quantity  varying 
linearly  with  time,  using  appropriate  scaling.   It  is  thus 

possible  to  study  the  kinetics  of  steady-flow  processes. 

2 
An  example  is  the  study  of  catalytic  isomerization  $  the 

results  of  mechanisms  assuming  either  single-site  or  dual 

site  adsorption  are  readily  obtained  for  a  wide  variety  of 

initial  conditions  and  flow  rates. 

A  computer  of  the  sort  described  is  of  only  moderate 

expense  -  a  quite  satisfactory  combination  need  cost  no  more 

than  a  good  polarimeter  -  and  can  be  operated  by  a  single 

chemist  or  chemical  engineer  working  in  his  own  laboratory. 

A  problem  of  the  type  described  can  be  set  up  and  solved  in 

a  few  hours. 
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